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I. INTRODUCTION 



The study of open charm production in high-energy lepton-hadron and hadron-hadron 
colhsions is considered as a test of the general apphcabihty of perturbative quantum chromo- 
dynamics (QCD) and also provides information on the parton distribution functions (PDFs) 
of protons and photons. The present analysis is to explore our potential to access a new 
dynamical regime, namely the high-energy Regge limit, which is characterized by the con- 
dition ^ /i ^ Aqcd, where is the total collision energy in the center-of-mass (CM) 
reference frame, Aqcd is the asymptotic scale parameter of QCD, and fi is the typical energy 
scale of the hard interaction. In the processes of heavy-quark production, one has /i > m, 
where m is the heavy-quark mass. In this high-energy regime, the contribution from partonic 
subprocesses involving t-channel parton (quark or gluon) exchanges to the production cross 
section can become dominant. Thus, the transverse momenta of the incoming partons and 
their off-shell properties can no longer be neglected, and we deal with Reggeized t-channel 
partons. 

The quasi-multi-Regge-kinematics (QMRK) approach is particularly appropriate 

for this kind of high-energy phenomenology. It is based on an effective quantum field theory 
implemented with the non-Abelian gauge-invariant action, as suggested a few years ago 
Our previous analyses of charmonium and bottomonium production at the Fermilab 
^ demonstrated the advantages of the high-energy factorization scheme over the 
collinear parton model as far as the description of experimental data is concerned. These 
observations were substantiated for 5-meson production at the Tevatron in Ref. js], where 
the experimental data were again well described using the Fadin-Lipatov effective Reggeon- 



Tevatron 



Reggeon-gluon vertex [2]. In Ref. 6|, where the effective photon- Reggeon-quark vertex was 
obtained and for the first time, the hypothesis of quark Reggeization was successfully used to 
describe experimental data on single prompt-photon production and on the proton structure 
functions F2 and Fl. 

The CDF Collaboration measured the differential cross sections da/dpT for the inclusive 
production of D^, D^, D*^, and Df mesons in pp collisions in run II at the Tevatron as 
functions of transverse momentum {px = \pt\) in the central rapidity (y) region j^]. These 
measurements were compared with theoretical predictions obtained at next-to-leading order 
(NLO) in the collinear parton model of QCD 8|, l9| taking into account quark and hadron 
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mass effects, and it was found that the latter improve the description of the data. 

The differential cross sections da/dpT arid da/dy for inclusive D*^ and Df photoproduc- 
tion measured by the HI 10(] and ZEUS 
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12| collaborations at the DESY HERA Collider 
were compared with NLO predictions in the collinear parton model. For D*^ mesons, this 
was done in three approaches: the zero-mass variable-flavor-number scheme (ZM-VFNS) 

nn n 

|13l. |14|. the fixed- flavor- number scheme (PENS) ^5|, and the general- mass variable-flavor- 
number scheme (GM-VENS) [l6|. The experimental results were found to generally lie above 
the NLO expectations. Eor Df mesons, the calculations were performed in the EENS |l^ 
and in the model suggested by Berezhnoy, Kiselev, and Likhoded |l7l |. 

In this paper, we study D-meson production under HERA and Tevatron experimental 
conditions as well as the charm structure function F2^c of the proton for the first time in 
the framework of the QMRK approach complemented with the quark-Reggeization 

hypothesis. This paper is organized as follows. In Sec. [Tll we present the basic formalism 
of our calculations and briefly recall the QMRK approach in connection with the quark- 
Reggeization hypothesis. In Sec. IIIIl we consider the charm structure function F2^c of the 
proton and compare our results with experimental data. In Sees. HV] and |Vl we describe 
Z^-meson production via c-quark fragmentation at HERA and the Tevatron, respectively. In 
Sec. IVIt we summarize our conclusions. 



II. BASIC FORMALISM 



In the phenomenology of the strong interactions at high energies, it is necessary to de- 
scribe the QCD evolution of the PDEs of the colliding particles (hadrons or photons) start- 
ing from some scale /io which controls a non-perturbative regime up to the typical scale 
fi of the hard-scattering processes, which is typically of the order of the transverse mass 
Mt = \J + p\ of the produced particle (or hadron jet) with (invariant) mass M and 
transverse momentum px- In the region of very high energies, which corresponds to the 
so-called Regge limit, the typical ratio x = /i/ becomes very small, x ^ 1. This leads to 
large logarithmic contributions of the type [a^ ln(l/a;)]"', where is the strong-coupling con- 
stant, which are conveniently resummed in the Balitsky-Eadin-Kuraev-Lipatov [18] formal- 
ism by the evolution of unintegrated gluon and quark distribution functions $^'^(x, /i^), 
where x and qx are the longitudinal-momentum fraction and transverse momentum of the 
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Reggeized parton w.r.t. the parent particle, respectively. Correspondingly, in the QMRK 
approach [s], the initial-state t-channel gluons and quarks are considered as Reggeons, 
or Reggeized gluons (R) and quarks (Q). They carry finite transverse momenta Qt with 
respect to the hadron or photon beam from which they stem and are off mass shell. 



The advantages o: 
factorization scheme 



the QMRK approach in comparison with the conventional krp- 



19| include: firstly, it uses gauge-invariant amplitudes and is based 



on a factorization hypothesis that is proven in the leading logarithmic approximation; sec- 
ondly, it carries over to non-leading orders in the strong-coupling constant, as recently proven 
{2^. The Reggeization of amplitudes provides the opportunity to efficiently take into ac- 
count large radiative corrections to processes in the Regge limit beyond what is included in 
the coUinear approximation, which is of great practical importance. 

Recently, the Feynman rules for the induced and some important effective vertices of the 
effective theory based on the non-Abelian gauge-invariant action [sl have been derived in 
Ref. 0- 

However, these rules only refer to processes with Reggeized gluons in the initial 
state. As for t-channel quark-exchange processes, such rules are still unknown, so that it is 
necessary to construct effective vertices involving Reggeized quarks using QMRK approach 
prescriptions in each application from first principles. Of course, a certain set of Reggeon- 
Reggeon- Particle effective vertices are known, for example for the transitions RR — > g |22| . 



QQ 
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23| . and RQ q |2J]. The effective 7*Q q vertex, which describes the 



production of a quark in the collision of a virtual photon with a Reggeized quark, has been 
recently obtained in Ref. 61]. 



In our numerical calcu 
Martin, Ryskin, and Watt 



ations below, we adopt the prescription proposed by Kimber, 



25| to obtain unintegrated gluon and quark distribution functions 



26|. 
26|. 



for the proton from the conventional integrated ones, as implemented in Watt's code 
To obtain the analogous unintegrated functions for the photon, we modify Watt's code 
As input for this procedure, we use the Martin-Roberts-Stirlin-Thorne 27|| proton and the 
Gliick-Reya-Vogt jasj photon PDFs. 



III. CHARM STRUCTURE FUNCTION F2,o OF THE PROTON 



On the experimental 
by HI y and ZEUS 



side, the charm structure function F2^c of the proton was measured 
SOj in deep inelastic scattering (DIS) of electrons and positrons on 
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protons at HERA. In this section, we consider this quantity in the framework of the QMRK 
approach endowed with the quark- Reggeization hypothesis. We thus need the partonic cross 
section for the production of a c quark in the colhsion of a virtual photon and a Reggeized 
charm quark. The relevant vertex was found in Ref. \^ and reads: 



iT - -r—^io + 



(1) 



where the four-momenta of the virtual photon, the proton, the Reggeized charm quark, 
and the outgoing charm quark are denoted as gi, P2, q2 = X2P2 + q2T, and k = qi + q2, 
respectively. We concentrate on photons with large virtuality = —qf ^ '^cj-.^^ that the 



massless approximation for describing DIS structure functions is appropriate 
obtain the following master formula for F2^c'- 



where Cc = 2/3 is the fractional electric charge of the c quark and X2 = xb{Q'^ + ^2)/<5^, 
with xb being the Bjorken variable. For definiteness, we choose the factorization scale to be 

In Fig. [H we compare the xb distributions of F2^c for various values of with the HI 



6|. We then 



291] and ZEUS [SOj data. We find good agreement for all values of Q^, except for the highest 



one, = 500 GeV^, where our prediction is about 50% below the data. This disagreement 
shows the importance of higher-order corrections at large values of Q^, which are beyond 
the scope of our present study. 



IV. D-MESON PRODUCTION AT HERA 



On the experimental side, ZEUS measured the pr distributions of [Uj and Df [1 



mesons with rapidity 3J] \y\ < 1.5 inclusively produced in photoproduction at HERA I, with 
proton energy Ep = 820 GeV and lepton energy Eg, = 27.5 GeV in the laboratory frame, 
in the ranges 2 < pj- < 12 GeV and 3 < < 12 GeV, respectively. In this section, we 
compare this data with our QMRK predictions. At leading order (LO), we need to consider 
only three 2 — > 1 partonic subprocesses, namely Cp7 c for direct photoproduction and 



c and RpCy — > c for resolved photoproduction, where the subscript indicates the 



mother particle. 
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Exploiting the factorization theorem, the pt distribution of direct photoproduction takes 
the form 

x$^(xi,ti,/i2)|M(C,7-c)|2, (3) 



where 



with kr being the transverse momentum of the produced c quark. We evaluate the quasi- 
real-photon flux /^/e in Weizsacker- Williams approximation using 



nl 1 '^max I 2-^2 



(5) 



•^7 Vmin \Vmm Vmax 

where a is Sommerfeld's fine-structure constant, mg is the electron mass, Q^in = T^lxtjil — 
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and Qmax is determined by the experimental setup, with Q^ax = 1 GeV^ in our case 
12| . As for the c ^ D fragmentation function (FF) Dc^d, we adopt the non-perturbative 



D*^ and Df sets determined in the ZM-VFNS with initial evolution scale /io = 3l| from 



fits to OPAL data from CERN LEPl. We choose the renormalization and initial- and final- 



state factorization scales to he fi = ymjy + p"^, where ttid is the Z)-meson mass. Using 
the Reggeized-quark-photon effective vertex from Ref . [6i] , the square of the hard-scattering 
amplitude is found to be 



|M(C7 ^c)|2 = 87rae24. (6) 

It is understood that also the contribution from the charge-conjugate partonic subprocess is 
to be included in Eq. 

In the case of resolved photoproduction via the partonic subprocess CpR^ — > c, the 
factorization formula reads: 

Pt-^ = J dy J dz J dx^ J dt2 J d(f)2fj/e{x^)z'^Dc^r){z,i/) 



X <l>?(xi,ti,/i2)$?;(x2,t2,/i')|A^(Cpi?^ ^ c)|2, (7) 



where 



Ptc^ ptc y 2 2 r Pt 

TT-FT, 2:2 = 7: —, ti =t2-2kT \/t2 cos <p2 + krp, t2 = q2T, = , 

Zzh/p Zx^zH/g z 



(8) 



with 02 being the angle enclosed between pt and q2T- Using the Reggeized-quark-Reggeized- 



gluon effective vertex from Ref. 



we have 



\M{CR^c)\'^ = -Tia,{l?)kl. (9) 

Again, the charge-conjugate partonic subprocess is to be included in Eq. ((Tj). Resolved 
photopro duct ion via the partonic subprocess RpC^ — > c is treated very similarly. 

In Figs.[2](a) and (b), our results for D*"^ and Df mesons, respectively, are broken down to 
;he Cj£y — > c, CpR-^ c, and RpCy c contributions and are compared with the ZEUS data 
111 . u\ . We find that the theoretical predictions are dominated by direct photoproduction 
and agree rather well with the experimental data over the whole pt range considered. 

V. i:>-MESON PRODUCTION AT THE TEVATRON 

CDF 0] measured the pt distributions of D", D^, D*^, and Df mesons with rapidity 
\y\ < 1 inclusively produced in hadroproduction in run 11 at the Tevatron, with -s/S* = 
1.96 TeV. To LO in the QMRK approach, the factorization formula for the CpRp c 
channel reads: 

da 



Pt'^j^ = J dy J dz J dti J d(f)i z'^Dc^d{z, jj^ 



X <l>P{xi,h,fi'm{x2,t2,fi')\M{CpRp^c)\\ (10) 



where \M{CpRp ^ c)|2 is given by Eq. ([9]), 

Xi = 7=, X2 = — 7=-, t2 = ti-2 — V^lCOS02 + ^- (11) 

zvS zvS z z^ 

The result for the RpCp c channel is similar and has to be included in Eq. (JTOi) together 
with those from the charge-conjugate partonic subprocesses. 

In Figs.[3](a)-(d), our results for D^, D*"^, and Df mesons, respectively, are compared 
with the CDF data [7]. We find that the theoretical predictions generally agree rather well 
with the experimental data, except perhaps for the slope. In fact, the predictions exhibit a 
slight tendency to undershoot the data at small values of px and to overshoot them at large 
values of pt- However, we have to bear in mind that these are just LO predictions, so that 
there is room for improvement by including higher orders. 
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In the framework of the coUinear parton model, comparisons with the experimental data 
of Ref. 0] were performed beyond LO, namely in the fixed-order-next-to-leading-logarithm 



(FONLL) scheme 



|32|. The FONLL predictions sys- 



and at NLO in the GM-VFNS U 

' ^ n 

tematically undershoot the CDF data [7]. The GM-VFNS predictions of Ref. [9|, which 
are evaluated with FFs determined in the ZM-VFNS [31j], describe these data within their 
errors, but are still somewhat on the low side. The degree of agreement is further improved 



32l | by evaluating the GM-VFNS predictions of Ref. 



9| using FFs extracted 



33| from a 



global fit to B- and Z-factory data of e^e annihilation in the very same scheme. 



VI. CONCLUSIONS 



In this paper, we explored the usefulness of the quark- Reggeizat ion hypothesis in the 
framework of the QMRK approach by studying several observables of inclusive charm pro- 
duction at LO, namely the charm structure function F2^c of the proton measured at HERA 



30l | as well as the one-particle-inclusive cross sections of D*^ and Df photoproduction in 



ep collisions at HERA ll|, [l^ and of D^, D^, D*^, and Df hadroproduction in pp collisions 



at the Tevatron [7]. In all three cases, we found satisfactory agreement between our default 
predictions and the experimental data, which is quite encouraging in view of the simplicity 
of our LO expressions for the partonic cross sections. By contrast, in the coUinear parton 
model of QCD, the inclusion of NLO corrections is necessary to achieve such a degree of 
agreement. We thus recover the notion that the QMRK approach is a powerful tool for the 
theoretical description of QCD processes in the high-energy limit and automatically accom- 
modates an iniportant class of corrections that lie beyond the reach of the coUinear parton 
model at LO U. 
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FIG. 1: F2^c{xb,Q^) as a function of xb at (a) Q 



(f) 500 GeV^. The HI Q (open circles) and ZEUS 



25, (b) 30, (c) 45, (d) 60, (e) 130, and 



30(] (filled circles) are compared with LO 



predictions from the QMRK approach with the quark-Reggeization hypothesis. 
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FIG. 2: pt distributions of inclusive (a) D*^ and (b) Df photoproduction for = 300 GeV and 
\y\ < 1.5. The ZEUS data from (a) Ref. [11] and (b) Ref. 3 

are compared with LO predictions 
from the QMRK approach with the quark-Reggeization hypothesis. 
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FIG. 3: pt distributions of inclusive (a) D^, (b) D^, (c) D*^, and (d) hadroproduction for 
\/5 = 1.96 TeV and \y\ < 1. The CDF data from Ref. [7] are compared with LO predictions from 
the QMRK approach with the quark-Reggeization hypothesis. 
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